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Abstract 

By considering the vacuum polarization, we study the effects of geometry on elec- 
trostatic self-energy of a test charge near the black hole horizon and also in regions 
with strong and weak curvature in static two dimensional curved backgrounds. We 
discuss the relation of ultraviolet behavior of the gauge field propagator and charge 
confinement. 
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1 Introduction 



One of the interesting areas in physics, is the study of behavior of classical and quantum fields 
in curved background and investigating how their properties are affected by the curvature 
of the space-time. It also provides a lot of insights into important problems such as black 
hole entropy, Hawking radiation, quantum theory of gravity and so on. 

It is well known that the change in the geometry of the space associated with the gravi- 
tational field, deforms electromagnetic field, inducing a self-force on a point charge at rest 
in a static curved space-time [1]. 

Also the presence of boundary condition like the boundary condition at conical singularity 
produced for example by a cosmic string or a point mass, alters the electromagnetic field of 
a point charge which after subtracting the infinite part, leads to a finite self-force [2]. 

A renewed investigation has been appeared in this subject in order to study the upper 
bound on the entropy of charged object by requiring the validity of thermodynamics of black- 
holes. This problem is studied in [3] , for classical black- hole backgrounds, in the absence of 
dynamical fermions, i.e. disregarding vacuum polarization. 

Another subject of studies in gauge field theory is the screening and confinement of 
charges. The static potential between external charges, which can be obtained from the 
Wilson loop expectation value, carries important information of infrared behavior of gauge 
fields which is suggested to be responsible for confinement, binding the quarks and anti- 
quarks into qq pairs (infrared slavery). Because of computational hurdles in four dimensions 
one can consider these problems in lower dimensional models, as a laboratory to study 
physical effects which can be carried out to the real world. 

In this letter we study the influence of curvature on the self-energy of static charges, 
by considering the effect of vacuum polarization in two dimensional static space-times. We 
show that in order to explain the confining behavior of QED2 [5] on a curved background, 
considering ultraviolet behavior of gauge fields and self-energy of external charges is neces- 
sary, in other words our method of studying confinement involves the behavior of two point 
function in the ultraviolet regime instead of the infrared. To do so We use the criterion 
expressing that in the confinement phase the energy of an isolated quark is infinite. In usual 
calculations in four dimensional flat space-time the self-energy of a test charge is infinite 
and is subtracted from the potential energy, but on curved space-times the finite part of 
self-energy is not a constant and must be considered in computing the forces [2]. In two 
dimensions this self-energy is related to the Green function of a Sturm- Liouville type op- 
erator at coincident limit and therefore is analytic. We show that, for QED2, the self-force 
can prohibit a single charge to be in some region of the space unless it is coupled to another 
opposite test charge, forming mesonic qq structure. We also obtain electrostatic self-force 
using the heat kernel method up to the adiabatic order four. 
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2 Geometrical effects on charge confinement and mesonic 
structure 



A general static two dimensional surface can be described by the metric 

ds 2 = y/g(x)(dt 2 - dx 2 ), (1) 

where y/g(x) is the conformal factor. 2 On this space-time, QED 2 consisting of charged mat- 
ter field interacting with an abelian gauge field in two dimensions is described by Lagrangian 

L = v^aWV (V, - ieAJi/;) + j^)^ ( 2 ) 

where 7 M are the curved space counterparts of Dirac gamma matrices. V M is the covariant 
derivative, including the spin connection, acting on fermionic fields, e is the charge of 
dynamical fermions. The dual field strength F, is described through F = e IIU d^A u , where 
^ u = ^iiv an d e 01 = — e 10 = 1. By integrating out matter fields one can obtain one loop 
effective action for the gauge field [4] 

where /i — In static case and using the Coulomb gauge A ± — this Lagrangian reduces 
to 

Hence in the presence of vacuum polarization, the gauge field has gained a mass via a 
peculiar two dimensional version of Higgs phenomenon. As a consequence, one may expect 
the replacement of the Coulomb force by a finite range force. We introduce two static 
opposite charges located at x = a and x = b, described by the covariantly conserved current 

J°( x ) = -s(x-a)), ^ = 0. (5) 

y9\ x ) 



The gauge field's equation of motion is 

d 1 dA 



H 2 A = e'(5(x-b) -8(x-a)). (6) 



dx y/g(x) dx 

The Green function of the elliptic operator ^ — ^ 2 satisfies 

, d 1 d 



" dx y/g(x) dx 



-H 2 )G{x,x') = 8(x,x'). (7) 



2 Euclidean version of this space-time is ds 2 — y/g(x)(dt 2 + dx 2 ) , whose one of the geodesies is the straight 
line parallel to the x axis. 
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In terms of G(x,x f ) the energy of external charges is obtained 



/2 

E = j T$dx = - J ' L eff dx = -^[G{a,a)+G{b,b)-2G{a,b)\. (8) 

This is the energy measured by an observer whose velocity = (g -1 / 4 ^), 0) is parallel to the 
direction of the time-like killing vector of the space-time and = ^-G(x,x) is the self-energy 
of a static point charge located at x. 

On flat surfaces this change of energy is 

E = —{l-e'^ b - a \). (9) 



This can be obtained by using the gauge fields Green's function 

G(x,x') = -— e-v^-A (10) 

Zi Li 

Note that —^-G(x,x), in contrast to higher dimensions has a finite value. This is due to 
finiteness of Helmholtz Green's function in coincident limit in one dimension. For distant 
charges as a result of screening, interaction term becomes zero, and E tends to self-energy 
of test charges, which in flat case is E^^^ — |- [6]. Note that although this self-energy 
is finite, but it is a constant therefore the self-force is zero. As we will show on a curved 
surface, the Green function of the Sturm-Liouville operator (7) is an analytic function at 
coincident limit and self-force becomes position dependent. 

Besides, the role of dynamical charge is affected by the presence of the curvature 

To see this we write the equation (7) as 

1 d ' d ^(x)G(x,x')-^G(x,x') = 5 -^l, (12) 



(13) 



y/9(x) dx y/g(x) dx •■• y/g(x) ' 

where G(x,x') = s /g(x)G(x,x'). This equation can be rewritten as 

(tk \ i\a 1 t d^/g{x) dG(x,x') 1 d 2 G(x,x') 5(x,x') 

(R(x) - „ )G(x, x ) + + Mx) dx2 = 

which is equivalent to 

{R(X) _ + _J_ ( A|£) )( ^_J_ G( ,^ )) + ^ ( ^) = SM . (14) 

Hence in the strong curvature limit ^> /i 2 , the Green function is approximately unaf- 

fected by dynamical fermions. In other words if we assume the same boundary condition for 
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the gauge field in the presence and absence of dynamical fermions, the vacuum polarization, 
in contrast to the flat case, doesn't change the energy and confining phase of the system. 

To elucidate this subject and to emphasize how the ultraviolet behavior of QED 2 on 
curved space-time is concerned in confinement of test charges let us give an example. Con- 
sider the following space-time in conformal coordinates 

ds 2 = dt2 - dx2 . x > o (15) 

x m v ' 

where m = 2 — 1/(^ + 1); k ^ 0, —1/2, this is one of the classical solutions of two dimensional 
scale invariant gravity [7] with a curvature singularity. Here we assume that this is only a 
classical background for QED 2 . The homogenous solutions G h (x) of the equation (7) satisfy 

± x ™±G h (x)-f, 2 G h (x) = 0. (16) 



For m 7^ 2, by defining z = x 1 ™ and Gh(x) = x ™ we obtain 
2 d 2 u du Am - l) 2 4u 2 z 2 . 

z ^ +z Tz- { w^w + w^w ,u=0 - (17) 



Hence 



G h (x) = x 2 m < 



^l^ilW Cl-mS^ 2 ) 



|Vy (2-m) 

x l^|(V(^F ;rl_T )- 



(18) 



Therefore the Green function is 



2 !_ m /m-if-^xt 2 )^™-i(^a;< 2 ) ifm>2 



m 



_! ,(^x< 2 )Ar, m -i,(=^-x> 2 ) ifm<2, 



where £>(<) is the bigger (smaller) of x, x' and J, K are modified Bessel functions. This 
Green function Satisfies Dirichlet boundary condition at x = and at x = oo. At the 
coincidence limit the Green function is 

G(x,a;) = -— ^— x 1 - m / |m -i,( r -^- r a; 1 -T )tf x ). (20) 

v ' ; \m-2\ l ~ |V |m-2| ; ' — |V |m-2| ; v ; 

In terms of the scalar curvature R(x) = -j=^, this relation becomes 

In the strong curvature limit, or in regions where (|i?(x)| 3> fJ. 2 ), by considering the asymp- 
totic behavior of Bessel functions this relation becomes 

fcGM = _|i+I| ( ^£>)*, (22) 
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which is /i independent as anticipated, k is defined after the equation (15). In regions where 
the curvature is weak |-R(x)| <C /i 2 , we have 



Note that for m — 0, G(x,x) = —j^, which is the same as the result on a flat uncurved 
surface [6]. 

Now we show the relation of electrostatic self-energy and confinement. G(x,x) in (21) is 
an analytic non-constant function for x > (which may become very large in some regions, 
signaling, as we will show, a confining situation), and therefore we expect that the charge 
feels an electrostatic self-force, affected by the geometry of the space-time as well as the 
boundary condition imposed on the gauge fields. 

We use the fact that in the confinement phase the energy of an isolated quark is infinite. 

Assume that m > 2; the energy needed to locate a single charge e' in the region (R(x) ^ 
0) <C fi 2 or x ~ (in the coordinate (15)), is 



which tends to infinity. In other words there is a great repulsive force on an external charge 
near x = prohibiting to have single charges in this region, or the energy of an isolated 
charge in this region is very large. The same procedure occurs in the region x ~ 0, for m < 2 
(when — > oo). So in these regions following the equation (8) only charges forming 
mesonic structure may survive. These two opposite charges must be near together in order 
to obtain a finite energy for the system. Our criterion for confinement is not based on the 
behavior of the energy in the infrared (where the geodesic distance of external charges tends 
to infinity). 

In QED^ the dielectric constant of vacuum is larger than unity as a result of the screening 
effect due to vacuum polarization. If instead one consider a case in which the dielectric 
constant of the vacuum vanishes then due to antiscreening effect the energy of the system 
becomes infinite unless we add another opposite test charge to the system then this fictitious 
system is confining. A similar phenomenon occurs in dual bag model. In that case quarks and 
anti-quarks must form mesonic structure in order to avoid divergences in static potential, 
i.e only color singlets have finite energy because the divergence term appears as a multiplier 
of total external charges [8] . In our model the role of dielectric constant in confining phase 
is played by the metric components (see equation (2)). To find some relations between 
geometry and permittivity see [9]. 

The repulsive forces on single charges besides the curvature of the space-time is related 
to the boundary conditions imposed on the gauge fields in defining the vacuum of the system, 




(23) 




(24) 
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for example in equation (26), although R is small but we have a great repulsive force. These 
repulsive forces may also be arisen on a surface with constant curvature (adS or dS space- 
times, obtained for example by taking m = 2 in the previous example) [10] . Also Maxwell field 
theory (disregarding vacuum polarization) in 2+1 dimensional conical space-times, despite 
the null curvature of the manifold exhibits a repulsive force on charges [2]. 

In the previous example we considered a space-time with a naked singularity, in this part 
we study the electrostatic self- energy on a black hole background. 

For Maxwell theory on a four dimensional Schwarzschild black-hole a test charge near 
the horizon is repelled by an image charge inside the horizon. In these cases one must 
subtract the infinite parts to obtain a renormalized Green function or the finite part of self- 
energy. On a two dimensional static space-time, in contrast to the Maxwell theory in three 
and four dimensions, as we noticed (after equation (10)) and will discuss later, electrostatic 
self-energy in QED2 is a well defined function. 

In Schwarzschild coordinate, we consider a non-extremal two dimensional static black 
hole described by the metric 

ds 2 = f(r)dt 2 - JJ^dr 2 - ( 25 ) 
At the horizon r = h, f(h) = 0. In this coordinate the equation (7) becomes 

(^ r )(^) -^) G M = /Wr,r'). (26) 

Near the (bifurcate) horizon, i.e. r ~ h, r > h, we have f(r) = n(r — h), where 2k denotes 
the surface-gravity. Assuming the gauge field tends to zero at infinity and is well behaved 
at the horizon, the two point function of the gauge field becomes 




G(r,r') = 2(r> - h)*{r < - h^K^J^ )h(M ~ )» ( 2 ?) 



where ?"<(>) is the smaller (bigger) of r and r' . 

To obtain the effect of the gravitational field on the electrostatic self-interaction, we use 
the global method used in [12]. If in a free falling coordinates the work 5W is needed to 
displace the charge slowly by a distance Sr, then this energy computed at asymptotic infinity, 
by considering the gravitational red-shift, will be (the space-time is flat at infinity) 



5E = sff{r)5W. 



Using the total mass variation law of Carter [11], and assuming that the metric is unperturbed 
by the presence of the charges [12], we arrive again to the equation (8). Note that the integral 
must be taken over Tq from the horizon to infinity. The self-energy of a test charge near 
the horizon is then 



E self .(x) = e /2 (r - h)/f 1 (2^L^)/ 1 ( 2 ^L_2). (29) 
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In contrast to the four dimensional case, there is an attractive force on the test charge near 
the horizon, which by considering the asymptotic behavior of Bessel functions is independent 
of vacuum polarization: this attractive force may be related to an image charge inside the 
horizon, when these charges are near together, that is near the horizon, the effect of vacuum 
polarization may be disregarded (self-energy is independent of fi). At the horizon the self- 
energy is zero, and in contrast to the previous example ultraviolet behavior of the Green 
function doesn't lead to charge confinement. Vanishing of self-energy may be understood 
as follows: Instead of the black-hole horizon one can consider an image charge inside the 
black-hole. Then the self-energy of the test charge is 

E = ^(G(r,r)-G(r,r')), (30) 

where r' is the location of the image charge — e', and G is the Green's function which does 
not satisfy the Dirichlet boundary condition. In the limit r = r' = h we obtain E — 0. 

Besides, if we assume r — > oo, the interaction energy between the charge and its image 
becomes zero, and we obtain E — |-, which is the self-energy of a test charge on a flat 
surface. This can be seen explicitly by setting f(r — > oo) = I. 



3 Heat kernel expansion of electrostatic self— energy 

In this part we study the short distance behavior of Green function G(x,x'), using the heat 
kernel of positive elliptic operator O := — + f^ 2 - This method can be used for slowly 

varying metrics. We write the heat kernel in the form 

h(r; x, x') = 7r= ex P( o fi 2 r)a n (x, x'), (31) 

n=0 V47T 2t 

which satisfies 

^, . ,. dh(r; x, x ) , . 

Oh(r; x, x') + V ' - ; = 0, (32) 
or 

where a — \ \ f*, g J (y)dy | 2 , is one half of the square of geodesic distance between (t, x) and 
(t, x') and r is the proper-time parameter. G(x, x') is given by 

[•CO 

G(x,x') = - h{T;x,x')dT, (33) 
Jo 

provided 

h(0;x,x') = 8(x,x'). (34) 

Therefore 

I oo _ 

G(x,x>) = --= Y,(^ +k) K nH (^)a n (x,x>) (35) 
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and for a = 0, 

1 OO 1 

G(x, x) = -— £ (^r n -^(n + -)a n (x, x), (36) 

2 V 7r n=0 2 

which is regular. 

Under a scale transformation parameterized by the positive number A, 

y/g{x) -> \y/g{x) 

R(x) - ^(x), (37) 

the Green function becomes Gx(/j, 2 ) = XG(/j, 2 X) or jG\(^) = G(jj, 2 ). We have written the 
\x dependence of G explicitly. Using (36 ), we obtain 

a„(A) = ?4^>. (38) 

Hence As a polynomial, a n consists only of rath power of g (including also its derivatives), 
where m — \ — |. In a n , the order of derivatives is In. For example as we will see, in 
ai(x,x), only the terms g~^(x)g' 2 (x) and g~%(x)g"(x) are present. 
In order to obtain heat-kernel coefficients we use the relations 

[a] = [a'] = 0, [<r< 2 >] = ^(x), [<7< s >] = ^=(x). (39) 

We have shown a(x, x) by [a] and / denotes the first derivative and (n) the nth derivative 
with respect to x. 

By solving the equation (32) for the Seeley coefficients, we obtain a recursion relation 

- -^9~^{x)g'{x)a' n {x,x') + g~^(x)a (2 \x,x') - (n + l)a n+1 (x,x') + 

-g'(x)a'(x,x')g~^(x)a n+1 (x,x') - g~^a'(x, x')a' n+1 (x, x') = 0. (40) 

For n < 0, a n<0 = 0. In order to satisfy (34) we must have a (x, x') = gs(x)gz(x'). Taking 
the diagonal value of (40) yields 

-\g-H*)sf (*)[<] +9-h*)[<>S ) ] ~ (« + l)K + i] = 0. (41) 

For n = 

N-^W/WKl + sifl, (42) 

hence 

N = -^<T Wto + 1 -g'Hx)g^(x). (43) 
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For [a 2 ] we require the diagonal part of ai derivatives: [a'J and [a^]. Differentiating (40) 
with respect to x gives 

( - \g-'hx)g^\x) + \g-ig'\x))W n \ - ^* (*)</(*) [a£ 2 >] + <T*[a<?>] - 

(n + 2)K +1 ] + ^- 1 (x)y'(x)[a n+1 ] = 0. (44) 

For n = 

K] = -\g-kx)g'(x)g^(x) + / + ^g-Hx)g^(x). (45) 

Another differentiation of (40) with respect to x in the limit x — > leads to the following 
equation 

( - ^r*O«O0 (3) (x) + V^y(^ (2) (*) - ^(s)</ s (*))k] + 



(-^ 2 \x)g-'Hx) + V(x)p"f (x))[ai 2 )] + </-*[<#>] " §</(s)<T* (*)[«?>] + 



2" v y " v y 4 
9 

4 

(^(x)^" 1 ^) - ^/(^- 2 (*))[a n+1 ] - ^(s^tok+i] - 

(n + 3)[a£Jj=0. (46) 

Therefore we find 

[ai 2) ] = -^9~Hx)g'(x)g^(x) + (x)# (4) (:r) - j^g~Hx)g^ 2 (x) - 

l ^g^(x)g'\x) + ||r^)/(^ (2) (*). (47) 

Equations (41), (45 ), (47) yield 

r i 245 15 , .2, v col / \ 2343 19 , .4 26 n , . . fji/ x 

N= 4 to WrW-^s 4 Ms w~i92^ 'WsWs 1 w - 

^-^s)^*) + Lg-i( X )gW( X ). (48) 

One can continue this method to obtain other [a n ]. Heat kernel coefficient can be expressed 
in terms of the scalar curvature R = ^g^^(x)g^(x) — ^g^^(x)g' 2 (x) and k(x) = \g'(x)g~ l (x) 
(at the horizon of a black-hole, 2k is the surface -gravity). 

M = \g^i.x)R(x) - ^g'^(x)K 2 (x) 

N = ^g-Hx)RW(x) - ±g-i(x) K (x)R!(x) - ^{x)R 2 {x) + 
^g-\( x )K 2 ( x )R(x) - ^-g-\{x) K \ x ). (49) 
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Hence the self-energy of the charge e' up to the fourth adiabatic order is 

e' 2 i 1 i 3 

Eseif.{x) = —g-^{x)[^ l g^{x) + -iT 3 (g*(x)R(x) - ^ 2 {x)) + 

A^(Ii2( 2 )(x) - k(x)R'(x) - \gHx)R\x) + 

f A n\x)R{x) -f A g-\x)n\x))]. (50) 

Note heat kernel expansion like the WKB method can not be applied at the horizon h or 
the turning point, i.e. where g(h) = 0. To obtain an expression for the self-energy near the 
horizon one can expand the metric and follows the steps after the equation (25). However 
the expansion (50), is consistent with the asymptotic behavior of electrostatic self-energy 
obtained near the horizon. For example Considering asymptotic behavior of Bessel function 
for large arguments in equation (29) gives the first term in equation (50). This first term is 
also the same as the result obtained for the small curvature limit (23). For y/g = 1, (50) is 
reduced to the flat case result E sel f. = x [6]. 
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